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Abstract 

We discuss a multistep variational approach for the study of many-body corre- 
lations. The approach is developed in a boson formalism (bosons representing 
particle-hole excitations) and based on an iterative sequence of diagonaliza- 
tions in subspaces of the full boson space. Purpose of these diagonalizations 
is that of searching for the best approximation of the ground state of the sys- 
tem. The procedure also leads us to define a set of excited states and, at the 
same time, of operators which generate these states as a result of their action 
on the ground state. We examine the cases in which these operators carry 
one-particle one- hole and up to two-particle two- hole excitations. We also ex- 
plore the possibility of associating bosons to Tamm-Dancoff excitations and 
of describing the spectrum in terms of only a selected group of these. Tests 
within an exactly solvable three-level model are provided. 

PACS numbers: 21.60.Jz, 21.10.Re, 21.60.Fw 
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1. INTRODUCTION 



Developing reliable microscopic approaches for the description of correlations in quantum 
many-body systems is a field of active research in various branches of physics. A preeminent 
role in this field has been traditionally played by the Random Phase Approximation (RPA) 
|T|]. Over the years, however, several attempts [^-[T^ have also been made to overcome the 
natural limitations of this theory related, in particular, to its lack of an internal consistency. 

In a recent publication |T^, with reference to the /3-decay physics and working in a 
quasi-particle formalism, we have discussed an approach aiming at improving the quality of 
the standard Quasi-particle RPA calculations usually made in this field. The basic point 
of this approach has been that of searching first for the best approximation of the ground 
state. In order to do this we have started with an initial Ansatz for this state and we have 
tried to improve this approximation through a series of minimizations which modified the 
structure of the state at each step. Excited states have, then, been constructed by acting 
with an RPA-like phonon on the ground state and fixing the variables of this phonon via 
the minimization of the energy. Tests of the procedure have proved to be quite encouraging. 
In the present work, however, we present an evolution of this method which we believe to 
be more effective and simpler to apply in realistic cases. 

The approach is developed in a boson formalism. As a preliminary step, then, a boson 
space will be defined where bosons identify particle-hole excitations and a mapping procedure 
will allow the transformation of fermion operators onto their images in this boson space. 



Similarly to the previous work ||T6|, the basic point of the approach will consist in searching 
first for the best approximation of the ground state. Differently from the mentioned case, 
however, this will be achieved by means of an iterative sequence of diagonalizations in 
subspaces of the full boson space. As a further and important difference from the case of 
Ref. |jl6[, as a result of this sequence of diagonalizations, a set of operators will also be 
generated which by acting on the ground state of the system will define a set of excited 
states. We will consider two cases: the case in which these operators carry only one-particle 
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one-hole (Ip-lh) excitations and the case in which they include up to two-particle two- hole 
(2p-2h) excitations. In both cases the comparison with exact calculations within a schematic 
model will allow us to judge the quality of the approximations. As a schematic model we 
have chosen the SU(3) model |r^ -pJ|]. 

In the second part of the work, we will reformulate the fermion-boson correspondence 
and identify bosons with properly chosen collective particle-hole excitations. Then, what is 
more interesting, we will explore the possibility of describing at least partially the spectrum 
of the system in terms of only a selected group of these bosons. Once again we will first 
search for the best approximation of the ground state and we will then construct excited 
states carrying up to Ip-lh and 2p-2h excitations. Tests within the schematic model will be 
provided also in this case. 

The paper is organized as follows. In Sect. 2, we will describe the basic points of the 
procedure limiting ourselves to the case of Ip-lh excitations. In Sect. 3, we will provide 
some applications within the SU(3) model. In Sect. 4, we will consider an extension of the 
procedure to include 2p-2h excitations. In Sect. 5, we will examine the case of bosons as 
collective particle-hole excitations. Finally, in Sect. 6, we will summarize the results and 
give some conclusions. The Appendix A will be reserved to describe some details of the 
mapping procedure. 



2. THE PROCEDURE: IP-IH EXCITATIONS 

To simplify the notation, we will illustrate the procedure directly within the exactly 
solvable model which has been used for our tests. This model, the so-called SU(3) model. 



was first discussed by Li et al. [0 and has been used more recently by Matsuo et al. 
and Takada et al. [|19| to test some approximation schemes. The model consists of three 
2r2-fold degenerate single-particle shells which are occupied by 2Q particles. In the absence 
of interaction, then, the lowest level is completely filled while the others are empty. This 
state, the "Hartree-Fock" (HF) state of the system, is denoted by |0). A single-particle state 



3 



is specified by a set of quantum numbers {j,m), wliere j stands for the shell (j=0,l,2) and 
m specifies the 2Q substates within the shell. The creation and annihilation operators of a 
fermion in a state (j, m) are defined by aj^ and ajm, respectively. 
Let us consider the operators 
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Kij = al^aj^ = 0, 1, 2). (1) 

m=l 

These operators satisfy the Lie algebra of the group SU(3) 

[Kij, Km] = SjkKii - SiiKkj. (2) 



It is assumed that the Hamiltonian of the model is written in terms of the generators K,, 



only and contains up to two-body interactions. Its form is |19 



i=l,2 «J=1,2 

+ J2 Vy{t,j,k){K,oKjk + KkjKo,), (3) 

i,j,k=l,2 

where the coefficients are real and obey the symmetry conditions V^(i,j) = Vx{j,i) and 
Vv{i,j) = Vv{j,i). The eigenstates of Hp are constructed by diagonalizing it in the space 

\n,n,) = -^={K,,nK,or\0) I , (4) 

where N'nin2 normalization factors. 



As in Ref. |(T^, we will work in a boson formalism. To begin, then, we define the boson 
space 

B=i\n,n,) = ^l=={b\r{blr\0)] , (5) 

where the operators b] obey the standard boson commutation relations 

[k,b]] = Sij, [h,bj]=0 (6) 
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and |0) is the boson vacuum. As evident from a glance at and (|^), a one-to-one cor- 
respondence exists between the states of F and B, the boson operators bj playing the role 
of the excitation operators KiQ and the boson vacuum |0) replacing the HF state |0). As 
anticipated in the Introduction, however, in Sect. 5 we will also examine a different corre- 
spondence and so a different meaning to attribute to these boson operators. 

The mapping procedure to construct boson images of fermion operators is the same 
discussed in previous works [Il2| , |l6i and it is based on the requirement that corresponding 
matrix elements in F and B be equal. The procedure is, therefore, of Marumori-type. We 
will give further details in App. A. Here, we simply say that, in correspondence with the 
Hamiltonian Hp (^, we introduce an hermitian boson Hamiltonian Hb which contains up 
to five-boson terms. This has therefore the general form 

HB = a + Y.(3^ibl + h.c.) + ^ %,blb, + ^ (f),,{blb] + h.c.) 

i ij i<j 

+ E E ^ijk{blb]bk + h.c.) + E E ^ijkiblb]bkbi 

i<j k i<j k<l 

+ E T.PiMblb]blbi + h.c.) + E Y.^ijkim{0Abihm + h.c.). (7) 

i'^j'^k I i<j<kl<m 

To illustrate the iterative sequence of diagonalizations on which our approach is based 
we start by introducing an arbitrary boson state I^E'o'''*)- We consider this as a zeroth-order 
approximation of the ground state and we assume |^(°)) = -i=(|0) + 6t|0) + 4|0)). Let's then 
consider the space 

5« = {|v^J°)),6t|v,(0)),6,|v,(0))^^ (8) 
and diagonalize Hb in this space. 

5(1) 

is, in general, considerably smaller than the full 
boson space B. In our calculations, for instance, we have assumed 2Q = 10 and this implies 
that the space B can have up to ten-boson states while B^^^ contains only up to two-boson 
states. However, if l'^'^^) denotes the lowest eigenstate resulting from this diagonalization, 
one can only expect that |\E'o^'') will provide an approximation of the ground state better 
than (or, at worst, equal to) l^t'o*''*). This is due to the fact that we are allowing the new 
state to have more components than I^E'o'''') and that the coefficients of I^E'o^'') are fixed to 
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guarantee the lowest energy of the state. We define I^I/q ) ^^e first-order approximation of 
the ground state. 

As a next step, let's consider the space 

(2) 

and diagonalize Hb in this space. If ) is the lowest eigenstate resulting from this 
diagonalization, the above arguments lead us to expect that also I^E'o ) ^ill be better than 
l^o^'')- ^6 define I'^'o^) the second-order approximation of the ground state. The procedure 
can go on as many times as one wishes. By performing a sequence of diagonalizations in 
spaces whose dimensionality remains unchanged (and much smaller than that of the full 
boson space) one can construct approximations of the ground state which improve step-by- 
step. 

It turns out to be interesting to reformulate the procedure just described as follows. Let's 
define the operator 

i i 

It is, then, 

\<^) = iQmo), (11) 

with ^ = = ^ (the coefficients remain undetermined in this case). Similarly, 
one can define an operator {QoY^^ such that 

i^s'^) = (QS)^'^i^r) (12) 

(k)\ 

and so on for all other approximations. In general, if |\l/o ) denotes the k-th approximation 
of the ground state, one can write 

l^f'^) = iQm<'-'^) = iQir\Qir''^ ■ ■ ■ iQino)- (13) 

Therefore I'^'q^) is a product of A; + 1 operators of the type (pUj), k corresponding to 
the k diagonalizations in the S*^*^) subspaces plus the operator (Qo)'-'^'* corresponding to the 
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starting Ansatz I^I/q ). Concerning this state, some comments are necessary to justify its 
use. In principle, one could have started with a diagonalization similar to all the other ones, 
namely in a space of the type where I^E'o^'') = However, the coefficients f3i of the 
boson Hamiltonian (|^) are nothing but the matrix elements of Hp between the HF state 
|0) and the Ip-lh states KiQ\0) (see App. A). These coefficients turn out to be zero in our 
model and the same would happen in a realistic case. In consequence of that no mixing is 
possible between the states |0) and 6]|0) and so a diagonalization in the space {|0),6]|0)} 
could generate (what indeed happens in our model) the boson vacuum |0) as the lowest 
eigenstate. This would lead to a crash of the iterative mechanism. 

Once a sufficient number of iterations has been performed the procedure is expected to 
reach convergence. If this is the case, any diagonalization beyond a given one, let's say the 
fc-th one, will have to leave the results unmodified. This necessarily implies that the operator 
{QoY''^^'^ which will emerge from the {k + l)-th diagonalization will have coefficients 

Convergence of the procedure therefore means convergence towards these values of the co- 
efficients X, Y and Z. 

As a result of the same {k + l)-th diagonalization, besides the operator {QqY'^'^^\ one 
will also obtain the operators (QI)'"'''^^^ associated to the remaining eigenstates. The number 
of these eigenstates is (up to) 2N, where N is the number of the Ip-lh excitations {N=2 in 
our model). If we call \gs) our best approximation for the ground state, i.e. \gs) = |\£'o'^'*), 
these further eigenstates can be written as 

= (^ = l,...,2Ar). (15) 

Keeping in mind that (Qq)^^'^^^ = ±1, one has that 

i9s\{QlY'^%s) = {z = l,...,2N). (16) 

Moreover, the operators (QI)*-'^"''^^ satisfy the orthogonality conditions 
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{9s\{Q^)^'^'\Q]f^'^\9s) = S,, (z,j = 1,...,2N). (17) 

This procedure therefore leads us to define a set of operators {QIY^'^^'' whose action on 
the ground state gives rise to a set of excited states which, considering the nature of the 
operators ([T0|), all carry excitations of the type Ip-lh. 

This way of representing the excited states shows evident similarities with that of RPA. 
The operators (p^ ) indeed remind us (although in a boson formalism) of the phonon 
operators of RPA. However, important differences do appear between the two approaches. 
While in RPA the ground state is defined as the vacuum of the operators Q, this is not true 
in this approach. Here, a sequence of operations of variational type leads us to construct 
both the ground state and the set of operators which define the excited states. In RPA, 
instead, the operators are first constructed by solving some equations (not of variational 
type) and an explicit expression for the ground state can be subsequently derived. We also 
notice that the Q^'s (p!0|) , although constructed in terms of bosons, are not real bosons 
themselves since they do not obey standard commutation relations of the type (^. This 
is not the case in RPA where, at least in the standard quasi-boson approximation, the 
operators are treated as bosons. 

Before concluding this section, some comments are necessary about the violation of 
the Pauli Principle which is always a risk whenever dealing with boson transformations. 
Although the sequence of diagonalizations described above can be extended as long as one 
wishes and so one can in principle form states of the type (jl^) which involve any number of 
operators Q\ not all the components of these states may be "physical" . In other words, there 
could be components, the so called "spurious" components, which have not a counterpart 
in the fermion space F. In our model, for instance, whenever acting with more than 2Q 
operators on the boson vacuum |0) one would form states having, among the others, 
components with more than 2Q operators 6] and these are all spurious components. 

In order to properly take into account these components, one should in principle perform 
the diagonalizations in spaces of the type 
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Bik) ^ {\^t\TBbl\^t\b,\¥t'^)}^^^^^, (18) 

where we have introduced the identity operator Ib of the boson space B. This is rather 
simple to do in the model under discussion. However, as we will see in the next sections, 
the rate of convergence of the procedure in the cases examined has always been such as to 
make the problem associated to the occurrence of these spurious components unimportant 
in these calculations. 

3. RESULTS 

The calculations we are going to describe refer to the following choice of the parameters: 
2Q=10, e(l) = e, e(2) = 1.5e, V,{i,j) = -2x, V,{i,j) = \x and Vy{i,j,k) = -fx (i, j,A;,= 
1, 2). Both e and x are parameters expressed in units of energy. 

In Fig. 1, the solid lines show the ground state energy (A) and the excitation energies of 
the lowest five states (B), in units of e, as functions of the strength x/e. These results are 
obtained by diagonalizing Hp in F. Dot-dashed lines show the equivalent results for Hb in 
B. The agreement between the fermion and boson spectra guarantees the very good quality 
of the boson image Hb- 

In Fig. 2, lower part, we show the ground state energies corresponding to different 
orders of approximation as indicated by the numbers which label the dot-dashed lines. For 
comparison, we plot (solid line) the energies which result from the diagonalization of Hb 
in B since this represents the best one can hope to reproduce in this approach. The same 
will be done in all the next figures. As seen in Fig. 1, however, fermion and boson energies 
differ very little from each other. A clear improvement of the quality of the approximation 
is observed in correspondence with the increasing of its order. 

Still in Fig. 2, upper part, dot-dashed lines show the spectrum obtained within this 
approach (the spectrum is found in correspondence with the best approximation of the 
ground state shown in the lower part of the figure). For comparison we also show the 
energies of the two one-phonon RPA states (dashed fines). RPA undergoes a collapse as 

9 



soon as the ground state energy starts deviating significantly from zero (x/e ~ 0.024). Tlie 
same states, but witliin tlie wliole range of x/e? are obtained within our approach and they 
well reproduce the exact ones. It is worth noticing that our approximated spectrum is 
formed by four excited states (with the only exception of very small values of the strength 
xl e where they can become two) as opposed to the two states of RPA. Already for a strength 
x/e ^ 0.014 the second one-phonon RPA state actually corresponds to the third excited state 
within the present approach. The forth approximate state lies higher in energy and has not 
been reported in the figure. 

4. 2P-2H EXCITATIONS 

The same procedure discussed in Sect. 2 and involving only Ip-lh excitations can be 
extended in a natural way to include 2p-2h excitations as well. The basic difference consists 
in performing each diagonalization of the iterative sequence in spaces of the type 

= {\^t\h\\^t\M\<"'^^^^^^^^ (19) 

Moreover, differently from the case of Sect. 2, there is no more need for an initial Ansatz 
for the ground state. The iterative procedure simply begins by performing a diagonalization 
m a space 5(1) of the form (|19D where \^f) = |0). 

To evaluate the role of these additional excitations in the structure of the spectrum, we 
have performed calculations similar to those shown in Fig. 2. The new results are reported 
in Fig. 3. For what concerns the ground state energy one observes a very good agreement 
and a faster convergence with respect to the Ip-lh case. Concerning the spectra of the lowest 
five excited states, besides the two "one-phonon" states discussed in the previous section, 
also the remaining states are now well reproduced (as for Fig. 2, this spectrum refers to 
the best approximation of the ground state as indicated in the lower part of the figure). As 
expected, then, the inclusion of 2p-2h excitations considerably improves the quality of the 
approximate spectrum. 



10 



5. BOSONS AS COLLECTIVE PARTILCLE-HOLE EXCITATIONS 



When performing the boson mapping we have estabhshed a one-to-one correspondence 
between the states \n1n2) and \n1n2) defined in the Eqs. (||) and (||), respectively. In such 
a correspondence, bosons 6] are images of the Ip-lh operators KjQ. However, as aheady 
anticipated, this is not the only possibility of correspondence. To show an alternative choice. 



let's proceed as in Takada et al. \T^\ and let's first define the Tamm-Dancoff (TD) phonon 
operator 



The amplitudes f,-^^ satisfy the orthogonality condition 

E-f^-r=^AV (21) 

i 

and are obtained by diagonalizing Hp in the basis {:^7=-ft'io|0)}j=i,2- In terms of these TD 
operators we construct the space 

This space is the same as (H) but just a different representation. Therefore, if we establish a 
one-to-one correspondence between states (^2|) and (|^) and we reconstruct the boson image 
of Hp, the new boson Hamiltonian will have different coefficients (for instance, the matrix 
7jj of (1^ will now be forced to be diagonal) but its spectrum will remain unchanged. In this 
new representation, bosons correspond to collective particle-hole excitations and so play a 
role very similar to that of the standard s,d,... bosons in the Interacting Boson Model picture 
pO| (where they are meant to represent collective particle-particle excitations). As in this 
case, then, it is natural to expect that the structure of the low-lying part of the spectrum 
may be described in terms of only a selected group of collective bosons. 

Our model appears particularly suited to illustrate this point. One can form two TD 
excitations (pO|) and their energies are shown in Fig. 4. As one sees, while increasing the 
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strength x/e, one of the energies remains almost constant while the other one shows a 
regular and sizeable decrease. This behaviour leads us to beheve that the lowest boson may 
play a preeminent role in the structure of the low-lying spectrum. We have made some 
calculations involving only this boson and (to start) only Ip-lh ecitations. By denoting 
6^(6) the creation(annihilation) boson operator associated to the lowest TD excitation, the 
procedure consists now in an iterated sequence of diagonalizations in spaces of the type 



As initial Ansatz we have assumed the state ) = ^(1^) + ^^l^))- '^^^ results are 
shown in Fig. 5. The agreement for the ground state and the first excited state is very 
good. It is worth stressing that this result has been obtained by working in spaces whose 
dimensionality is remarkably smaller than that of the full boson space B (3 versus 66 in our 
model) . 

In Fig. 6, we show similar calculations which involve up to 2p-2h excitations. The quality 
of the agreement for the ground state and the first excited state (already good) remains 
basically unchanged while also the lowest "two-phonon" state (to use an RPA language) 
is now well described. These calculations therefore confirm the expectation that only the 
boson associated to the lowest TD excitation plays an active role in the structure of these 
states. 



In this paper we have presented a multistep variational approach for the study of many- 
body correlations. The approach has been developed in a boson formalism (bosons repre- 
senting particle-hole excitations) and based on an iterative sequence of diagonalizations in 
subspaces of the full boson space. Purpose of these diagonalizations has been that of search- 
ing for the best approximation of the ground state of the system. The procedure has also 
led us to define a set of excited states and, at the same time, of operators generating these 




(23) 



SUMMARY AND CONCLUSIONS 
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states as a result of their action on the ground state. We have considered two cases: a), the 
case in which these operators carried only Ip-lh excitations and, b), the case in which also 
2p-2h excitations were included. 

The approach has been tested within an exactly solvable three-level model. The compar- 
ison betwen exact and approximate ground state energies has allowed us to appreciate the 
convergence of the procedure in both cases a) and b). In the first case, a comparison also 
with standard RPA calculations has shown that the "one-phonon" states that this theory 
could reproduce up to its crash-point were now well reproduced in the whole range of vari- 
ation of the strength. With the inclusion of 2p-2h excitations also the remaining low-lying 
states of the spectrum have been well reproduced. 

In the second part of the paper, we have reformulated the fermion-boson correspondence 
and identified bosons with Tamm-Dancoff phonons. We have then explored the possibility 
of describing at least partially the spectrum of the system in terms of only a selected group 
of these bosons. In our model this has implied restricting the set of two possible bosons 
to the one corresponding to the lowest excitation energy. We have verified that the ground 
state was still well reproduced and so was the first excited state already at the level of Ip-lh 
excitations. 

The possibility of selecting a restricted set of collective particle-hole excitations and 
therefore of constructing the boson space only in terms of the corresponding bosons appears 
to be quite appealing. It may represent, in fact, an effective way to reduce the dimension- 
alities of the system and so to lead to a much simplified application of the procedure to 
reahstic cases. 
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APPENDIX A: 



In Sect. 2, we have established a one-to-one correspondence between a set of fermion 
states F = {\n1n2)} (Q) and a set of boson states B = {|rii?22)} (J^), both sets being orthonor- 
mal. In correspondence with a given fermion operator Op-, the mapping procedure which we 
adopt searches for a boson operator Ob such that matrix elements between corresponding 
states be equal. The operator Ob defines the image oi Op va. B. 

The construction of proceeds step-by-step involving at each step matrix elements 



between states which belong to increasingly larger subspaces of F and B ||T2|. Having in 
mind to construct boson images with no more than five-boson terms, it is sufficient to involve 
at most the two subspaces 



F' = \0),^=Ki,\0),^=K,oK,o\0),^==K,oK^oK,o\0) 



B' = l\o)A\o),^bmi^i=bm\o)\ , 



and 



where J\ff, Af^^, Nfji^, and A/^^^ are normalization factors. The boson image Ob which 
one derives is an hermitian operator which has the form (|^ and coefficients: 

a= (0|6i.|0), 
(0|dFif.o|) 



Ml 



_ (O|iro,0^ir,o|) 



^ (O|0^ir,oi^,o|) 



Sijkl — 



where 



and 



_ {Q\KoiKo,KokOFKio\0) _ ^ , . (3) 



'<j' 



_ {0\KoiKojKoi.d,KioKrno\0) 



B 
Im 



Z^i' Pi'^iji,i'lm + Y,i'<j' ^i'j'l^iji^i'j'm + ^i'<j' ^i'j'm^ijk,i'j'l 



^ijk,i'j'k' — i^ii'^jj'^kk' + ^ii'^jk'^kj' + ^ij'^jk'^ki' + ^ij'^ji'^kk' 

+Sik'Sji'6kj' + ^ik'^jj'^ki') /J^^k- 
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FIGURES 

FIG. 1. Ground state energy (A) and excitation energies of the lowest five states (B) as functions 
of the strength x/e- Solid lines are obtained by diagonalizing Hp (^) in F (^) while dot-dashed 
lines refer to Hb @) va. B 

FIG. 2. Ground state energy (A) and excitation energies of the lowest five states (B) as 
functions of the strength x/e- Solid lines are obtained by diagonalizing Hb @ in i? (^) while 
dot-dashed lines are obtained with the procedure described in Sect. 2 (only Ip-lh excitations). 
Numbers label different orders of approximation. Dashed lines (B) show the RPA one-phonon 
energies. 

FIG. 3. The same as in Fig. 2 but dot-dashed lines refer now to calculations involving up to 
2p-2h excitations. 

FIG. 4. Tamm-Dancoff energies as functions of the strength x/e. 

FIG. 5. Ground state energy (A) and excitation energies of the lowest five states (B) as 
functions of the strength x/e- Solid lines are obtained by diagonalizing Hb @ in B (^) while 
dot-dashed lines refer to calculations involving only the lowest Tamm-Dancoff boson and Ip-lh 
excitations. Numbers label different orders of approximation. Further details in Sect. 5. 

FIG. 6. The same as in Fig. 5 but dot-dashed lines refer now to calculations involving up to 
2p-2h excitations. 
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Fig. 1 
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Fig. 2 




Fig. 3 
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Fig. 5 




Fig. 6 
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